Constraining the reionization history with large 
angle CMB polarization 



Loris P L Colombo 

Dipartimento di Fisica "G. Occhialini", Universita di Milano-Bicocca, Piazza della 

Scienza 3, 120126 Milano, Italy 

I.N.F.N., Via Celoria 16, 120133 Milano, Italy 

E-mail: loris . Colombo @mib . inf n . it 

Abstract. The first-year WMAP data release showed that the reionization optical 
depth to CMB photons is greater than previously thought. This follows from 
unexpectedly high values of the Cf E spectrum, for I up to ~ 30, presumably allowing 
a measurement of the E-mode polarization spectrum Cj , sooner than expected. This 
note aims to test the capability of large-angle polarization experiments to explore the 
history of the cosmological reionization, considering also the impact of cosmic variance. 
In particular, here we discuss how well the ionized fraction x e and the reionization 
redshift z r can be separately measured, at various levels of instrumental noise. 
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1. Introduction 

In the last ten years, measures on Cosmic Microwave Background (CMB) have reached 
high sensitivity and resolution (see, e.g., [1, 2, 3, 4, 5]), providing the spectrum of 
temperature anisotropies, Cf, up to / ~ 1000. This allowed to better constrain various 
parameters of the cosmological model. Among them, the optical depth r, due to 
reionization, has a peculiar character. Other parameters directly reflect initial conditions 
and/or inflationary outputs; the value of r, instead, conveys information on the physical 
conditions in the Universe due to more recent events, and opens an observational window 
on the physics of ancient cosmic objects. Furthermore, while other parameters influence 
data through linear physics, the cosmic opacity r is clearly linked to later non-linear 
evolution. 

Unfortunately, r is only mildly constrained by Cf, because of a degeneracy with 
the primordial spectral index (n) of density fluctuations ([6, 7]; see also [8]). For this 
reason, pre- WMAP Cf data could only constrain t to be ^ 0.40, at the 2-a level [9]. 
On the other hand, several models of reionization hinted at a value of the optical depth 
in the range 0.03 — 0.05 (see, e.g, [10]). 

The recent WMAP first-year release t has significantly modified the context, thanks 
to data on the Z?T-correlation spectrum, Cf T . Quite in general, Cf T and the E- 
polarization spectrum Cf contain independent information on r. In particular, the 
most part of this information is conveyed by the spectral components with / ^ 30; it is 
their unexpectedly high amplitude that led to estimate r ~ 0.17 ± 0.04 [11]. 

In the presence of such signal, it is legitimate to wonder whether current (or 
more advanced) experiments can succeed in providing us with more information on 
reionization, besides the value of r. In this work we shall refer to the expected 
performance of the Sky Polarization Observatory (SPOrt:): [12, 13, 14]), first to outline 
its capacity to provide an independent confirm of the high r value, by using high- 
sensitivity polarimeters operating at three MW frequencies (22, 32 and 90 GHz). Then, 
we shall focus on the capability of low-/ polarization experiments to inspect the cosmic 
reionization history, by analyzing the expected outputs of experiments run at increasing 
levels of sensitivity. Our conclusions directly apply to any (nearly-) full sky experiment 
with a design similar to SPOrt; using them for possible future experiment, with minor 
technical differences, is substantially harmless. 

At large angular scales, cosmic variance (CV) can never be neglected and risks to 
reduce the strength of CMB data in constraining the reionization history. One of the 
main difficulties of this analysis amounts to taking carefully into account its effects. We 
shall see that, as these effects interfere with low-/ inspections, a sensitivity much higher 
than the one of WMAP or SPOrt ought to be reached, before the full potentiality of 
large angles results is exploited. 

f http : / /lambda . gsf c . nasa . gov/product/map/m_overview . html 
| http://sport.bo.iasf.cnr.it/ 
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2. Monte Carlo analysis 

2.1. Model properties 

An optical depth r = 0.17 ± 0.04 corresponds to a reionization redshift z r ~ 17 ± 5, 
assuming a sharp reionization in a ACDM model with WMAP best fit parameters. 
On the other hand, the Gunn-Peterson effect seems to affect high-z quasar radiation 
([15, 16]), emitted at redshift ~ 6; this shows that at least a fraction (possibly as 
small as 1 %) of neutral hydrogen was still present at that time. Taken together, these 
observations imply that reionization is a fairly complex process, and that its complete 
characterization requires more than the value of the optical depth ([17]). 

To account for such an observational picture, several models have been proposed, 
including double reionization models ([18, 19]) or a period of extended reionization 
([20]), but more complex pictures are possible. A feature common to most these models 
is a bump in Cf and Cf E spectra at low /, but other effects can be expected, involving 
spectral components up to I > 1000. The analysis carried on in this work is based on 
a partial reionization model, in which baryonic materials reionize abruptly at a redshift 
z r , attaining a reionization fraction x e < 1. Thus, a given value of r corresponds 
to different x e — z r pairs. Although this is admittedly a toy model, it may reasonably 
approximate a large family of possible reionization histories. We shall take other relevant 
cosmological parameters as known from independent data (e.g., from high-/ spectral 
analysis); we fix their values as follows: a flat spatial section with reduced density 
parameters Q m h 2 = 0.147, Qbh 2 = 0.024 (for matter and baryons, respectively), the 
Hubble parameter is 70 km/s/Mpc and the primeval spectral index is n — 1. We assume 
no contributions from tensor modes of gravity perturbation; this means, in particular, 
no 5-mode polarization. In a previous work ([8]) we tested that, in practice, likelihood 
results are insensitive to reasonable variations of these parameters. 

Parameter estimation from CMB data usually follows a Bayesian maximum 
likelihood approach (see, e.g., [21, 22]). In this work, however, we are mainly interested 
in polarization measures at large angles, to test how they can improve the reconstruction 
of the reionization history. Therefore, we adopted a Monte Carlo approach to directly 
account for CV, and performed a wide set of realizations of four models, corresponding 
to r = 0.17 or 0.13 with x e = 0.6 or 1. 

Artificial anisotropy and polarization data were built using CMBFAST* and 
HEALPix + codes. The resulting maps have been smoothed with a purely Gaussian 
filter of 7° FWHM, to represent SPOrt instruments beam, and we chose a HEALPix 
resolution N side = 16, corresponding to a pixels' width of ~ 3.5°, thus approaching the 
Nyquist frequency of the configuration. At the angular scales considered here, the non- 
ideality of SPOrt's feed horns give rise to a spurious polarization signal due mainly 
to the temperature anisotropy signal smoothed on scales equal to the instruments' 

* http : //physics . nyu. edu/matiasz/CMBFAST/cmbf ast . html 
+ http : //www . eso . org/science/healpix/ 
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FWHM. Assuming a rms T intensity of ~ 30/jK, this translates into a rms contamination 
< 0.2/iK, which is roughly an order of magnitude smaller than the expected noise of the 
SPOrt experiment [12, 23]. 

We considered increasing levels of sensitivity for polarization (noise-per-pixel a P i X = 
2, 1, 0.5, 0.1 /xK). The first value should be achieved by combining the two 90 GHz 
SPOrt-channels, after two years of data taking, while using the 22 and 32 GHz channels 
to map and remove the galactic foregrounds contamination. Temperature noise was 
held fixed to o pix = 1/iK, roughly the value expected for the 4-year WMAP results. 
As on large scales signal-to-noise ratio in temperature signal is much more than unity, 
varying temperature noise in the same range of cr pix does not affect our results. Noise 
was assumed to be uncorrelated, both between different pixels and between the different 
Stokes parameters. For each choiche of r, x e and o~ P i X , 1000 independent sets of artificial 
data were considered. 

2.2. Likelihood evaluation 

When multi-parameter analysis for small angle experiments is required, likelihood 
computation rapidly becomes time intensive, forcing the adoption of various 
approximations (see, e.g., [22, 24]) or even a fisher-matrix approach. In this work, 
however, we are interested mainly in how polarization measures at large angles can 
improve the reconstruction of the reionization history. Moreover, SPOrt large beam 
size implies that data will be binned to a fairly low number of pixels. Together, these 
considerations make the evaluation of the full likelihood function in coordinate, or pixel, 
space feasible for a reduced number of parameters. This has the added advantage of 
simplifying the joint likelihood between experiments with different sky coverages. 

An artificial CMB sky map is a realization of a given cosmological model M. whose 
statistical properties (i.e. the sets of Cj coefficients, with Y = T, E, TE) are known. 
On each pixel, measures of temperature fluctuations (T) and of the Q and U Stokes 
parameters are ideally performed. 

Galactic foregrounds contamination is a serious issue in T maps, we therefore 
exclude the Galactic plane (i.e. the region of sky with Galactic declination \5 g \ < 20°) 
from anisotropy maps. In the case of polarization, Galactic contamination, while still 
an issue, plays a less important role. The first year WMAP data showed that the most 
prominent foreground contribution is due to Galactic synchrotron, whose polarization 
spectrum has a steep dependency on frequency. According to the synchrotron template 
developed by the SPOrt team [25], at 90-100 GHz CMB polarization should dominate 
over foregrounds contamination by a factor of 2-3 to almost 1 order of magnitude, 
depending on the value of the optical depth. Accordingly, we exclude from polarization 
maps only the equatorial polar caps, with declination \S\ > 51.6°, which SPOrt will 
not be able to inspect from its position on board the ISS. Notice that the polar caps 
partly overlap with the Galactic plane. The number of pixels available for anisotropy 
(polarization) measures is therefore N T = 1984 (N P = 2404). These data form 
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a (N T + 2iV P )-components vector x = {T(i = 1, N T ), Q(i = 1, N P ),U(i = 
1,....,N P )}. 



The statistical properties of the process underlining the synthetic data vector can 
be expressed through the correlation matrix (x T jXj) = = Sjj + N^-, where the 
brackets mean ensemble average. In particular, the correlation matrix has been written 
as the sum of a signal term S^, which depends on cosmology and geometry, and a term 
which accounts for detector noise. 

The expected anisotropy correlation between two pixels with angular separation $ij 
is given by 



where (cos $;.,•) are Legendre polynomials and Bf are reduction coefficients due to 
pixelization and beam smoothing. For other correlations, similar expressions hold and 
can be evaluated using the relations given in [21]. 

We assume no correlation among noise in different modes and pixels, so Ny is 
diagonal and has distinct values (o~p ix ) 2 and a^ ix in the former Nt and in the latter 2Np 
pixels, respectively. 

Assuming Gaussian statistics, the likelihood of a model M' given the synthetic 
data x reads then: 



In a Bayesian analysis, the posterior probability density is given by the product of the 
likelihood function and the prior probability density. According to Bayes theorem, in the 
case of an uniform prior the most likely model is the one for whom the likelihood function 
is maximum. Confidence regions in parameter space are then found by a likelihood ratio 
criterion. These regions are also called credibility regions and generally represent only 
an approximation to the true confidence intervals, particularly when CV is significant. 

A frequentist analysis, instead, stems from the question: What is the probability 
of obtaining a given output from an observation if the underlying physical process has 
a certain set of parameters? To answer this question, the usual procedure is to generate 
a large set of synthetic observations, analyze them as real data, and then study the 
distribution of the results. Confidence regions, which in this case are exclusion regions, 
are determined as the portion of parameter space which contain a set percentage of the 
total results. 

3. Results 

For each simulated sky, likelihood distributions have been evaluated on the x e -z r plane. 
We show results in two ways: first, we average over the whole set of 1000 realizations 
and show 1 and 2-<r confidence regions (see figure 1 and figure 2). Outside of them, we 
have exclusion regions, i.e. the portions of the x e -z r . that can be excluded with a given 
degree of confidence, taking into account the variance both in the detector noise and in 
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the true CMB data. For the range of a pix we tested, with the exception of the lowest 
value, the detector noise variance however dominates polarization measures, while the 
proper CV is more relevant for anisotropy data, where the signal to noise ratio is higher. 

However, when CV is significant, as it is for the angular separations considered 
here, the statistical properties of a realization can differ significantly from those of 
the underlying process. It is therefore significant to show how maximum likelihoods 
parameters are distributed. In figure 3, the distribution on r is shown, both for a pix = 2 
and 0.5 /iK; solid (dashed) lines refer to the case x e — 1 (= 0.6). This plot is meant 
to test how far r determinations are affected by the actual reionization history. Then, 
in figure 4, we report the maximum likelihood distribution over x e , after marginalizing 
over z r . Finally, in figure 5, the distribution over x e is shown, after marginalizing over 
z r . The same values of a pix are considered, while solid and dashed lines refer to the 
same cases as above. 

4. Discussion 

Figure 1 shows the 1 and 2-er confidence regions for models with r = 0.17, when 
op = 2^K. They clearly show a likelihood distribution elongated along the constant 
optical depth directions; this confirms that CMB data, for current noise levels, can 
provide constraints only on the integrated value r (see, e.g. [26], [27]), but also that 
such constraints are fairly safe. 

Figure 2 shows the 1 and 2-cr confidence regions, still for models with r = 0.17, 
but when a pix = 0.5/iK. At this pixel noise level, the polarization sensitivity allows us 
to constrain an additional parameter (x e or z r ) beside r, at least with a 68% c.l. Once 
again, by marginalizing, e.g., over x e , we obtain the behaviour of the likelihood as a 
function of r. Figure 3 confirms that determination of the optical depth is substantially 
independent from the reionization history, at least within the range of reionization 
models considered here (however, other reionization histories can bias r determinations 
[28]). This is true both for a pix = 2 and 0.5 (J.K. 

By inspecting figure 4, we notice that, for the noise levels currently achievable 
(a pix ~ 2 yuK), the distributions of maximum likelihood values are greatly superimposed 
and it is not possible to distinguish between the two reionization histories. On the 
contrary, for o P i X = 0.5 //K, the superposition is restricted to the extreme tails of 
histograms; accordingly, at such a noise level, some details on the reionization history 
can be inspected. Notice that these plots provide the probability at which a given model 
can be excludedby the available data, instead of telling us the likelihood of a given model 
being the correct one, as is instead done in Bayesian analysis. 

For both sensitivities, however, the distribution corresponding to x e — 1 is more 
sharply defined than the one corresponding to x e = 0.6. This is due to the fact that 
in the former case the true model necessarily lies on the boundary of the portion of 
parameter space considered. Analogous considerations could be made if, instead of 
marginalizing on the reionization fraction, we had marginalized over z r , as it is shown 
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in figure 5. 

As outlined above, in this work we inspected a number of further cases, for which no 
plot is given. In particular, models with r = 0.13 show quite a similar behavior, although 
some difference, on the level of noise required to achieve a clean separation of the two 
models, exists. In particular, the overlap between maximum likelihood histograms, for 
o~ P ix = 0.5 fiK worsens from 2% to 4.5 %. For all cases, considering a pix = 1 fiK yields 
intermediate results. 

We report here also results obtained for a pix = 0.1 fjK. While this noise level 
is certainly beyond the capabilities of current instruments, it should be achieved by 
the PLANCK§ experiment. While much of the information on reionization histories 
considered here is carried by the low-multipoles, PLANCK's higher spatial and frequency 
resolutions provide an advantage in the form of reduced sidelobes contamination and 
a better monitoring of foregrounds. Simulations show that clean separation between 
the two different reionization histories should be achieved at more than 95% c.L, even 
taking into account CV (see figure 6). Moreover, at these noise levels, contribution 
from S-mode polarization should no longer be ignored, and could provide additional 
constraints on reionization. These, too, should be detected by the PLANCK satellite. 

5. Conclusions 

After the WMAP first-year results, it has become clear that reionization history is 
a complex process and CMB measures can provide unique information about it. In 
particular, CMB polarization (and, to a lesser extent, temperature-polarization cross 
correlations) encode a great deal of information not only on the total optical depth, 
but also on the evolution with time of the ionized fraction. However, this knowledge 
is found mainly in the first I ^ 30 multipoles, i.e. at large angles. These scales are 
strongly affected by cosmic variance; therefore, the statistical properties of the actual 
sky can differ significantly from those of the underlying cosmological model. 

In this work we investigated whether measures of CMB power spectra at large 
angular separations can provide actual constraints on reionization history. In particular, 
we considered the SPOrt project as a benchmark of a large angles polarization 
measurements, though we also explored instrumental sensitivities better than those 
of the actual experiment. We considered toy models with the same optical depth, 
but different fractions of ionized baryons; in both case we assumed an instantaneous 
reionization. 

Using a frequentist analysis, we confirmed that, at the noise levels currently 
achievable, CMB data can provide informations mainly on the total optical depth. The 
SPOrt experiment will therefore provide the first independent confirmation of WMAP 
results on the optical depth, in a manner completely free of any spurious leakage between 
the T and E signals. 

§ http : / / astro . estec . esa.nl/Planck 
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Distinguishing between a model in which all baryons reionized from a model in 
which the ionized fraction reaches only 60% requires, instead, noise levels significantly 
lower. In particular we found that, for r = 0.17, discriminating between the two models 
at 68% c.l. requires a polarization pixel noise a p i x = 0.5//K, while for <: 95% c.l. 
o~ P ix ~ 0.1 fiK is needed. While the models of reionization considered here are admittedly 
simplified, they provide two extreme examples of behaviour while being consistent with 
the currently available data. Moreover, our main aim was the determination of the 
sensitivity needed to extract information on an additional reionization parameter from 
CMB data. A complete assessment of the effects of a generic reionization history on 
CMB power spectra is certainly beyond the scope of this note, as it requires a completely 
different approach. In particular the non-uniform, or patchy, nature of the reionization 
process should be considered. This usually implies additional features at multipoles 
/ 1000, which need to be inspected in order to fully exploit the information on 
reionization encoded in CMB. 

Nonetheless, our results confirm that large angles polarization data can be used 
to discriminate between different reionization histories, but that next generation CMB 
(polarization) experiments will be indispensable for shedding light on those details of 
the reionization process that can be inspected through this observational window. 
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Figure 1. Expected likelihood contours in the x e — z r plane for models with r = 0.17 
and two different reionization histories; a pix — 2/iK. Solid and short-dashed lines trace 
the boundary of the 68% and 95% confidence level regions, respectively. In red we 
plot results for a model in which reionization is complete, while blue lines refer to a 
model in which x e — 0.6. Likelihood contours are elongated along direction of constant 
optical depth end do not allow to distinguish between the two models. 
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Figure 2. Same as figure 1, but with a p t x = 0.5/xK. Likelihood function are elongated 
along lines of constant r, it is possible to distinguish between the two models at least 
at 68% c.l. 
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Figure 3. Likelihood function after marginalizing over x e , for the models of figure 1. 
Results are plotted as a function of r instead of z r . Solid lines refer to complete 
reionization, dashed lines to x e = 0.6. No significant bias in the determination of the 
optical depth is introduced by the reionization histories considered here. 




Figure 4. Histogram of the distribution of the maximum likelihood value of z r , 
for models with r = 0.17 and different reionization histories. Solid lines refer to 
x e = 1, z r ~ 17, dashed lines to x e — 0.6, z r ~ 24. Neat separation of the two 
distributions requires pixel noises ~ 0.5/xK. 
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Figure 5. Same as figure 4, but for x e . 




Figure 6. Same as figure 1, but with a P i X = 0.1/iK. Clean separation between the 
two models is achieved at better than 95% c.l. 



